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Abstract 

An object’s wake in a plasma with small Debye length that drifts across the mag¬ 
netic field is subject to electrostatic electron instabilities. Such situations include, for 
example, the moon in the solar wind wake and probes in magnetized laboratory plas¬ 
mas. The instability drive mechanism can equivalently be considered drift down the 
potential-energy gradient or drift up the density-gradient. The gradients arise because 
the plasma wake has a region of depressed density and electrostatic potential into 
which ions are attracted along the field. The non-linear consequences of the instability 
are analysed in this paper. At physical ratios of electron to ion mass, neither linear 
nor quasilinear treatment can explain the observation of large-amplitude perturbations 
that disrupt the ion streams well before they become ion-ion unstable. We show here, 
however, that electron holes, once formed, continue to grow, driven by the drift mech¬ 
anism, and if they remain in the wake may reach a maximum non-linearly stable size, 
beyond which their uncontrolled growth disrupts the ions. The hole growth calcula¬ 
tions provide a quantitative prediction of hole profile and size evolution. Hole growth 
appears to explain the observations of recent particle-in-cell simulations. 


1 Introduction 

The wake behind an object in a plasma that drifts perpendicular to the applied magnetic held 
is hlled in by plasma how along the held from either side. This how produces a characteristic 
multidimensional potential well structure in the wake that attracts ions and repels electrons. 
The approximate steady-state form of supersonic wake potentials of separated ion streams 
has been established through one-dimensional models for decades [H El [3] and more recent 
work has established the subsonic solution requiring multiple dimensions|H |5]. However, 
there remains a great deal of uncertainty about the stability of the wake to unsteady short 
wavelength electrostatic perturbations. The solar-wind wake of the moon [HI [71 El E] is a classic 
naturally-occuring example of this wake problem, and in-situ satellite measurements have 
observed various electrostatic huctuations in it mm- Several large-scale computational 
simulations [T^ IT3l [Til [15] have also shown wake instabilities, but their nature has been 
controversial. The purpose of the present work is to provide a detailed explanation of the 
mechanisms that drive the wake instabilities. These may have important applications also 
for magnetized laboratory plasmas and their interactions with probes. 
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The idealized configuration we study is represented by a magnetized plasma fiowing 
perpendicular to the field but normal to a fiat, thin, object [I6]. This is equivalent to a plasma 
fiowing with a sufficiently high Mach number past a spherical (or similar approximately unity 
aspect ratio) object. The high cross-field velocity in this second case, causes the object to 
be thin relative to the characteristic lengths in the wake. In other words, the sphere is 
strongly compressed in the flow direction, when measured in appropriately scaled units. 
The analysis represents the plasma velocity distribution function in one dimension, along 
the assumed uniform magnetic field. The plasma is presumed to drift in the transverse, 
wake direction, with simply a uniform drift velocity. So there is a one to one correspondence 
between downstream position and time since passing the object’s position. The Debye length 
is much smaller than the object. 

A self-consistent wake potential develops that attracts ions and repells electrons illus¬ 
trated in Fig. In a previous paper[T6] it has been shown by integration along orbits that 



Figure 1: Potential energy structure of electrons as a function of (x) distance down the wake 
and {y) position along the magnetic field. (Potential is cut off for viewing purposes at height 
5). Example collisionless orbits that arrive at a position on the wake axis are shown. Those 
that have low parallel energy at the final postion are substantially de-energized. Artificial 
electron/ion mass ratio of 1/20 is used to assist the visualization. 

the collisionless electron distribution in the wake potential structure acquires a depression 
that is localized in velocity. We here call this localized reduction of fe{v) the “dimple”. 
It is generated on electron orbits that are near the threshold of being reflected by the po¬ 
tential energy hill. These orbits climb the hill, converting their parallel kinetic energy into 
potential energy. Then, because they approach the peak with very small parallel velocity 
(they are nearly or just reflected) they spend a long time near the ridge of the potential, and 
during that time their transverse drift carries them down the potential ridge. Eventually 
their parallel motion carries them down off the ridge, but not before they have substantially 
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reduced their total energy compared with when they climbed it. They have experienced 
“drift de-energization”. By contrast electron orbits that are far from the threshold of re¬ 
flection, having either much more or much less total energy than the potential ridge, spend 
much less time on the potential ridge. They are far less de-energized. The distribution 
function is constant along orbits in a collisionless plasma. So if the external distribution is 
monotonically decreasing in kinetic energy (e.g. a Maxwellian) then an orbit that started 
(outside the potential structure) at a higher total energy (because of de-energization) has a 
phase-space density f{v) lower than orbits that have experienced less de-energization. This 
is the qualitative explanation of the mechanism forming the dimple. Its form was calculated 
quantitatively by numerical orbit integration in reference |16j . 

The distribution function dimple that arises is linearly unstable to Langmuir waves. 
Therefore one expects this de-energization effect to excite electrostatic instabilities, which 
will have a tendency to £11 in and smooth out the distribution non-linearly until the growth 
rate is suppressed. Because the dimple size depends strongly on the electron to ion mass 
ratio, the free energy available prior to non-linear saturation (which was calculated) also 
depends on mass ratio; simulations that use artificially low mass ratio are therefore liable 
to obtain unphysically large fluctuation levels. For true mass ratios the energy available 
to instabilities is only 10“^ to 10“^ of the electron thermal kinetic energy; so the level of 
Langmuir wave turbulence expected is modest. 

The purpose of the present work is to pursue further the non-linear development of the 
instability driven by this de-energization mechanism so as to explain what is observed in 
recent large-scale simulations of this problem|17j. Those simulations clearly observe the 
formation of the electron distribution dimple, but the observations are of course of its self- 
consistent non-linear state. In phase space (see Fig. [^, the dimple is centered along the 
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Figure 2: Illustrative electron phase-space orbits in the vicinity of the wake’s potential energy 
hill. The hill extends along the parallel y coordinate to the edge of the wake. The dimple 
lies along the lower left to upper right branches of the separatrix, hatched. 

separatrix contour of constant total energy that the wake potential defines. In addition to 


3 






























incoherent noise, we observe small persistent coherent structures, like eddies, localized to the 
dimple in phase space. They propagate in along the dimple with approximately the local 
electron velocity (and acceleration). These structures, which (we will show) are electron 
holes, therefore leave the simulation, moving to large parallel distance, quite quickly. The 
exception to this behavior is that holes move much more slowly near the x-point of the energy 
contours (where the phase-space velocity is zero) which is naturally at the ridge of the wake 
potential structure. As the simulations progress (down the wake) eventually one (or more) 
of the electron holes near the x-point grows to a large size, and disruption of the ion velocity 
distribution occurs. There is a large amount of free energy in the ion distribution, because the 
ions are in two streams, of modest energy spread, with opposite velocities (> c*). The strong 
hole growth and disruption of these streams occurs at a place where (time when) linearized 
calculation indicates that the ion streams are stable because of their large separation. The 
puzzle that our current analysis addresses is how the perturbation becomes large enough to 
disrupt the ion streams well before they themselves become linearly unstable. Our answer 
is that the mechanism is a non-linear one involving electron holes. 

In section we formulate the Boltzmann equation for the parallel electron distribu¬ 
tion, and solve it approximately analytically in a potential of specihed shape to hnd the 
electron dimple at the potential ridge. This solution supplements the prior numerical orbit 
integrations [16] by providing an analytic form for the dimple, in particular its velocity-width. 
In view of the substantial approximations required to achieve this analytic solution, section]^ 
approaches the problem instead, by an integration in the parallel direction rather than along 
the two-dimensional orbits. This alternative (and equivalent) formulation shows that the 
drift (de-energization) effects can be conceptualized as a term in the one-dimensional Boltz¬ 
mann equation of approximately the “Krook” collision form. Solving this equation gives an 
identical expression for the dimple, through a conceptually different set of approximations. 

The second formulation is more useful for incorporating the effects of presumed quasi- 
linear diffusion hlling in the dimple. Section explains the expected consequences of a 
self-consistent level of incoherent turbulence. It is shown that this system cannot explain 
the growth of perturbations to sufficient amplitude to disrupt the ions and tap into their 
energy until a place on the wake is reached where the ion streams are very close to linear 
instability. In other words, it cannot explain what is observed in the simulations. 

Section 1^ provides an explanation and analysis of the coherent electron holes, and shows 
that the drift de-energization mechanism can be equivalently regarded as the drift convection 
of holes into regions of higher background density. This effect causes holes to grow in depth 
and velocity width. The self-consistent growth of hole width with background density and 
the resulting hole prohle (for quasi-neutral holes) is calculated analytically for Maxwellian 
background electrons and beam ions. Holes that retain their integrity and remain near the 
ridge of the wake potential structure (the x-point) can grow to sizes sufficient to disrupt 
the ion streams when the density increase is of the order of one e-folding. Moreover their 
characteristics are consistent with what is observed in the simulations. Therefore we interpret 
the pre-linear-threshold disruption of the ion streams as caused by the long-term non-linear 
growth of electron holes until they become energetic enough to tap the ion free energy. 
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2 Solving For the Dimple in Two Dimensions 

2.1 Boltzmann’s Equation with drift and quasilinear diffusion 

Including ad hoc quasilinear velocity-space diffusion [18] with coefficient D, we can write 
Boltzmann’s equation as 



v.V-fa.vJ/ = V„(Z1V,/). 


( 1 ) 


which for the one-dimensional distribution, magnetized case with coordinate y in the mag¬ 
netic field direction, and x in the perpendicular direction becomes 


dt y cit ^ dy m dy dv J dx dv \ dv J ' 


( 2 ) 


Here the right hand side can be considered the additional source terms in the 1-D Boltzmann 
equation arising respectively from drift de-energization and quasilinear velocity space diffu¬ 
sion. We will consider a time-independent situation: ^ = 0 and constant drift Vx- Velocity 
V written without a subscript refers here to Vy, the parallel velocity, and the distribution 
function is one-dimensional along y. 

The orbits are the characteristics of the left-hand side. They are the paths in (parallel) 
phase space corresponding to constant energy 


S 


1 

2 


+ q(f) = const. 


(3) 


They are most easily found as the contours of constant energy in phase space. 


2.2 Collisionless Orbits in Specified Potential 

We take parameters to be normalized so that velocities are in units of the cold ion sound 
speed Cs = \jTe/mi, and potential is in units T^je. The perpendicular distance is scaled 
such that X =distance/M_i_ where M_i_ = Vxjcs- Then the normalized energy equation for 
electrons becomes E = \mrV^ — 0, where = mg/mj. 

The dimensionless potential form is considered to be controlled by dynamics separate 
from what happens to the instabilities. The specihc form illustrated in Fig. [Tjis based on an 
approximate solution in the form of two expansions of plasma into a vacuum, patched at the 
symmetry axis as discussed previously (TB]. However, the only features of this potential shape 
that substantially matter in the present context are that —0 is symmetric and single-peaked 
in I/, having a known curvature near the ridge at ?/ = 0 and being zero beyond a certain y- 
distance; and that it decays from a large value at a: ~ 0 monotonically in the perpendicular, 
i.e. downstream wake (x) direction. So we will simply specify that 

• 00 (^) is the monotonic potential at ?/ = 0; 

• the curvature is given in terms of a scale length w by = —(j)o/w‘^{x)-, 

• and y = Y[x) is the edge of the perturbed potential. 
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[A model potential that fits simulation results with T, = Te has 0o = —l/1.2a:, y = 1 + 2x, 
and w = (1 + 2a;)/1.66.] 

Now since in normalized units 


1 dy dy 
Cs dt dx' 


(4) 


the equation of the electron orbits in 2-D space may be written in the vicinity of the ridge, 
in terms of an expansion as 

d^y d(j) 00 ..X 

dx^ oy 

Because of the smallness of the orbits do not have a large duration (a:-extent). So it 
makes sense to approximate (po/w'^ as a constant, ignoring its x-dependence. Then the orbit 
can be solved trivially as 

Hr. 

( 6 ) 


Vo 


where 


y = — sinh(fc[x — xq]) 

K 


k'^ = —(po/rrirw'^ 


(7) 


and Xq and Vq are the position and velocity of the orbit when y = 0. The approximations 
leading to this expression are not well justihed near the edge of the perturbed region (and not 
at all outside it), nevertheless most of the orbit of electrons that cross the ridge slowly is spent 
near the ridge. It is therefore reasonable to use eq. (§ to estimate the orbit duration X = 
Xo — X (considered the duration either in time or in space, since x and t are interchangable) 
from the edge of the perturbed region to ?/ = 0. It can then be considered the solution of 


sinh(A;X) = kY/\vo\. 


( 8 ) 


In order for the approximations adopted to be consistent, both sides of this equation must 
be large compared with unity. Therefore sinh 4 gxp and 


kX ^ In 


2kY 

l^^ol 


In 


/m^w\vo\ 


(9) 


2.3 Resulting Distribution Without Diffusion 


The electron distribution function ai y = 0, when D = 0, can be deduced by considering the 
change in total parallel energy S arising from the perpendicular drift term. To the extent 
that y X applies to the relevant parts of the orbit (which we’ve already assumed to be a 
good approximation), the potential energy change arising from cross-field drift (which is the 
drift de-energization) can be estimated from the x-gradient of the potential at the ridge 

(to hrst order in X/xq). Substituting for X we get 


6£ ^ 


<700 

kdx 


ln(2fcy/|xo|) 


<700 

dx 



In 


/ 2/^y \ 

yy/m;w\vo\ j ■ 


(11) 
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This energy is not included in the parallel energy conservation along the orbit. In other 
words, denoting the kinetic energy at the start of the orbit by /Coo = and when it 

reaches the ridge /Co = \mrvl, we have 

/Coo = /Co ~ 00 + (12) 

Consequently ai x = xq, y = 0, the distribution function /o is different from the (presumed) 
Maxwellian at the orbit start: 


exp(0o) exp(-^mrno) 


exp (00 -lCo-5£) 

n P 

mrW\Vo\ 


2C|0o|V' 


(13) 


where 


d0o 

kdx 


d(j)o w , _ 

-- ^=^mr 

\/W 


(14) 


At the characteristic distance down the wake xq ~ 1, all quantities 0o, d(j)o/dx, w, and Y 
are of order unity. So P is of order the square root of the mass ratio Ynp. 



Figure 3: The dimple in the electron distribution at y = 0, /o as a function of normalized 
electron velocity u, for three different electron to ion mass ratios The cross marks the 
dimple width, where /g = 0. The parameters are 0o = —1.32, d(j)o/dx = 0.96, w = 1.15 (so 
P = 0.96ydTv), and Y = 2. 

The dimple in the electron velocity distribution function is illustrated in Fig. in whi ch 
electron velocity has been expressed normalized to its thermal velocity u = v^j \J 2 TeJnie = 

VQJmr/2. The dimple arises in this Vlasov case as a hnal multiplication factor on the 
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otherwise Maxwell-Boltzmann behavior. It actually dominates the behavior near uq = 0 
forming a cusp. Strictly speaking our approximations are quantitatively unjustihed as u —)■ 0, 
but the qualitative observation that a cusp forms is correct. The approximations also assume 
mrUo(= 1 and so the factor in brackets is small. But the power P to which the factor 

is raised is also small. The derivative of /o with respect to Vq from this expression can readily 
be shown to be zero where 

= -rurvl = -P (15) 

which is a small quantity, of order the square root of the electron/ion mass ratio, m^. This 
condition may also be written 


^^0 


d(f)o w 


dx 


\(po 


1/2 


1 / 4 ’ 

rrir 


(16) 


which indicates the approximate width of the dimple in velocity space. This width is dictated 
by the electron/ion mass ratio. Therefore simulations that use artihcially increased mass 
ratio (rrir = me/rrii 3> 1/1836) will increasingly misrepresent the electron behavior. 


3 Solving for / by Parallel Integration 


Thus far we have approached the problem accounting fully for two space dimensions and 
discussing integration along orbits in 2-D space plus 1-D velocity. Solving the problem 
analytically has required major approximations, but has given a reasonable estimate of the 
result when there is no velocity-space diffusion. 

A different approach to solving for the dimple is to do integration not along the 2-D 
spatial orbits but along only 1-D (y) in space. That is actually how eq. ([^ is organized. 
2-D orbit integration takes the hrst term on the RHS to be part of the orbit characteristics 
(and so far has not included the second diffusive term). By contrast 1-D integration of the 
equation leaves the U 3 ;-convective term on the RHS, and integrates along a hxed x path in 
which y and v vary. This is then a truly 1-D (but phase-space) treatment, but instead of / 
being constant on orbits, it varies in accordance with the terms remaining on the RHS. 

Within this perspective, we can regard the two terms on the RHS as being dimple¬ 
generating convective de-energization, and quasilinear diffusion. To some degree they will 
balance one another: one tending to form the dimple, the other to smooth it away. The 
solution for / at some position y can be found in principle by starting in the unperturbed 
background region yy (actually at // = R, the edge of the potential perturbation), and 
integrating orbits inward to position y. In principle the solution is simply 

f(v, V,) - nyy, = + dL (17) 


This integration must be taken along phase-space orbits of constant £ = -|- q(j){y), and 

in practice needs to be done in terms of position. 


f{y,Vy) - f{yY,VY) 



dl 

dx 


+ 


dv 



dy 

V 


(18) 





This integral determines the difference between the actual / and the Maxwell-Boltzmann 
approximation (a Maxwellian scaled by exp((;/))). 

A solution by this technique requires us to know what the value of the terms in the RHS 
integral are. Focussing hrst on the convective de-energization term Vxdf /dx, we don’t know 
its value exactly until we have the solution everywhere. However, it may in some circum¬ 
stances be reasonable to approximate it in a manner which avoids us having to solve the 
full-scale integro-differential system. One such approximation is to presume that the shape 
of the dimple changes only slowly with x-position. If so, then the dominant contribution to 
df /dx can be estimated to be the variation of the overall level of /, which is approximately 
the Maxwell-Boltzmann. Its variation with x (in normalized parameters) is / oc exp(0), in 
which case 


= -Vx. 


(19) 


Vx^ _ _ 
f dx dx 

Here is like a collision frequency. And indeed, the term in the Boltzmann equation to 
which this approximation corresponds is a “Krook” collisional term —Vxf- Since the orbits 
of interest spend most of their time near the potential ridge at ?/ = 0, it is reasonable to take 
Vx to be uniform, given by taking 0 = 0o. That is the main approximation of this treatment. 

Incorporating just this term (i.e. taking H = 0) for now, we can perform the integral 
along constant-x based upon the resulting equation 


1/ d ld(j)d\ 1 df 

f \ dy nirdydv) f dt 


( 20 ) 


whose solution is 

fiy,Vy)/fiyY,VY) = exp[-z/^(t - ty)] (21) 

Thus the deviation from Maxwell-Boltzmann distribution can be considered to be an expo¬ 
nential multiplicative factor whose argument is proportional to the time an orbit takes to 
reach the position y. 


Notice that the dimensionless form for d/dt in eq. (20) shows the terms on the LHS 
are usually very large compared with the convective de-energization term (the term). 
That term is important only where the LHS terms are nearly zero, i.e. near y = 0 where 
= 0 and at values of v nearly equal to zero. In other words, the dimple generation takes 
place predominantly at the axis, for velocities near zero there. However, it is not that the 
contribution to df /dt is larger there, it is that the orbit spends far more time there than 
anywhere else. Passing that region in phase space contributes most strongly to t — ty. 

The duration, X = t — ty, of the orbit to the position y = has already been solved for in 
It was there taken as an approximation that variation of df/dy with x could be 

Instead 


2.2 


section 

ignored. Here it is no approximation, because we are integrating along x =constant 


the approximation has been made in df /dx. In any case, we can immediately appropriate 
the solution eqs. ([8| and as 


X 


1 

k 


In 


2kY 
l^'ol , 



w In 


2^1001 H 

Jm/w\vQ\ 


( 22 ) 


It should be no surprise that substituting this result into eq. (21) gives exactly the same 


dimple as previously: eq. (13) 
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What we’ve demonstrated, therefore, is that the dimple formation can be calculated by 
explicit integration at fixed x, based upon an approximation of the drift de-energization 
term in a Krook form. This demonstration gives some additional conhdence in the prior 
treatment. But it also gives us a more direct way to incorporate quasilinear diffusion, and 
(later) to understand electron hole growth. 


4 Quasilinear Electron Velocity Diffusion 

Now we consider the effects of instabilities that will arise as the dimple forms and prevent 
it from ever becoming the deep cusp that the Vlasov treatment finds. One possible result 
of such instabilities, if they consist of many incoherent modes, is an effective quasilinear 
velocity-space diffusion. That’s the case we discuss first. 


4.1 Self-consistent Diffusion level 

Under this assumption, the physics of the steady state is that the diffusion magnitude, D, 
adjusts itself corresponding to a moderate time-independent level of turbulence sufficient to 
maintain the distribution function at an approximately neutral stability. It is reasonable (if 
the Debye length is small compared to other lengths in the problem) to assume that the 
growth rate of the electron instabilities is intrinsically large compared with other timescales 
in the problem. In that case, D must adjust itself so that the instability threshold is never 
signihcantly exceeded. In other words, marginal stability is always approximately satisfied. 

Solving for the distribution function in those circumstances requires us to suppose that 
as we integrate along the orbit we encounter levels of quasilinear diffusion that are just 
sufficient to maintain the distribution function marginally stable. Doing so requires that in 
phase-space regions where the RHS terms are important (i.e. mostly near y = 0 and u = 0) 
the diffusion term counterbalances the convective de-energization term. 

When the two RHS terms exactly balance 


When D is (approximately) independent of v, and / deviates only a small amount from con¬ 
stant, i.e. in the vicinity of a shallow dimple, the solution is of this equation is a parabola. 
However, we don’t require that the terms exactly balance for orbits whose duration is suf- 
hciently short that the perturbation introduced by is smallj^ The dimple can therefore 

^ One way to model that fact is to allow the product of the RHS times the orbit duration to be no larger 
than some appropriate quantity. For example, if we require no more than a modest fractional reduction P 
of / in the dimple, then we must take 


X{v) 



dv 



<Pf- 


( 24 ) 


Where X is large, the Pf term is small and we recover the previous condition. But for larger r;o, when P/X 
becomes comparable with , the Pf term decouples the diffusion term from any necessity to balance the Vx 
term, and it can subside to zero. All this seems rather more elaborate than justified by the current precision. 
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be considered to be constrained to have a positive second derivative (equal to Vx/D) over 
a region around n = 0 that extends to a speed (|n|) at which VxX becomes smaller than of 
order unity. Outside that velocity region, the second derivative of / can become negative, as 
it must in order to merge the dimple with the bulk of the electron distribution function. The 


details of that outer region depend on how quickly the constraint eq. (23) is relaxed, whether 
D varies with v and so on. Such details cannot be precisely calculated using the analytic 
principles on which this treatment is based. Some ansatz must be adopted. A simple and 
plausible one is to choose to represent the dimple as a negative Gaussian perturbation to 
the bulk Maxwellian distribution. The velocity width of the dimple Gaussian, expressed as 
2vd such that the Gaussian is oc exp(—is determined by UxX ~ 1. 

The plasma fluctuation level adjusts the diffusion coefficient D to achieve marginal sta¬ 
bility. If D is small, the dimple is deep, because its magnitude is such as to give second 
derivative VxjD near its peak. If D is large the dimple is shallow. Thus, orbit duration de¬ 
termines the width, and marginal stability determines the depth of the dimple. The dimple 
width is given (see eq. by Ud ~ vo ~ (in units of Cg) or Ud ~ (in units of 

Vte)- 


4.2 Electron Marginal Stability 

The dispersion relation of electrostatic waves ise = l-l-x = 0. Instability requires the real 
part of the susceptibility Xr to be negative at frequency uj in the upper half of the complex 
plane where the imaginary part of x is zero. A bulk Maxwellian electron distribution (ig¬ 
noring ions for now), contributes a susceptibility real-part approximately 3fJ(Xe) = ^/k'^^'oe 
(at wave phase velocities small compared with the electron thermal speed). The contribu¬ 
tion from a dimple Gaussian of temperature Td and negative density —Ud is the same but 
multiplied by —ndTe/rieTd- So the total (real part) electron susceptibility is 


^iXe) 


1 


ndTV 

neTd- ' 


(25) 


Because k is essentially a free choice, it can be adjusted for any negative value of 3fJ(Xe) to 
make 3?(Xe) = ~1- For a symmetric distribution such as we are considering, the imaginary 
part of the susceptibility, A(xi), is zero at a; = 0; so negative 3fJ(xe) is sufficient (as well as 
necessary) for instability. Marginal stability of electrons alone is therefore at ndTe/n^Td = 1, 
which means 


Ue 


Td 

= — ~ m 


1/2 


( 26 ) 


However, if we just focus on the zero-velocity peak. 


/d(0)_nd /Td\F^ 

/e(0) ne \Td) \tJ 


(27) 


The dimple depth at marginal stability should be a quite noticeable decrease in the distri¬ 
bution function at zero velocity, fractionally 1/1836^/^ ~ 1/7. Fig. [^illustrates some cases. 
Gomparison with Fig. [^ shows that these depths, if anything, somewhat overestimate what 
is expected from filling in the collisionless dimple by diffusion. 
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Figure 4: Electron distributions representing a dimple of negative Maxwellian form with 
temperature equal to that are marginally stable with immobile ions. 
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A crnde npper bonnd on the magnitnde of the potential pertnrbation that will produce 
quasilinear diffusion snfficient to maintain marginal stability can be estimated as follows. 
Qnasilinear diffusion presnmes the cnmulative effect of stochastic orbits produced by multiple 
modes of different phase velocities. How many modes are involved is uncertain, but what 
is certain is that it is at least greater than one, and that therefore the amplitude of any 
one unstable mode is insnfficient of itself to flatten the distribntion (at marginal stability). 
A mode whose phase-space island size is eqnal to the width of the dimple (~ 2vd) is large 
enough on its own to flatten the distribntion. The single-mode pertnrbation snfficient to 
create snch an island width is (in normalized nnits) (j) ~ and is the npper bonnd of 
the pertnrbed potential at qnasi-linear marginal stability. This potential makes a very small 
pertnrbation to the ions, a fractional energy pertnrbation of only m\l‘^ jvi < ~ 1/50. 

More generally, quasilinear velocity-space diffnsivity of particles by a resonant spectrum 
of waves of specihed electric held (or potential) is proportional to the inverse sqnare of the 
particle mass. The ion streams’ velocities place them inside the dimple, snbject to the same 
resonant spectrnm of waves as the electrons. They will experience a diffusivity smaller by 
a factor = 1836“^. Negligible ion pertnrbation occurs at electron marginal qnasilinear 
tnrbulence levels. The free energy of the ions cannot be tapped by qnasilinear electron 
instabilities. 

4.3 Ion susceptibility contribution 

When the ions have effectively a two-stream distribntion in the region under consideration, 
they contribnte further to instability by negative contribntion to the real part of the snscep- 
tibility. The contribntion for low-temperatnre eqnal-density beams of velocity ±Uj (in nnits 
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of Cg) is 

This is sufficient to make a system with purely Maxwellian electrons unstable when n* < 1, 
which is the upper edge of the ion-ion instability region in a “Stringer” plot jl90 If the ion 
speed is not a great deal higher than this threshold, then the ion contribution modihes the 
electron marginal stability condition, rendering it: 

^ = ^(1 - l/v!) ~ my^(l - l/vf). (29) 

The marginal-stability depth of the dimple is somewhat decreased. And to make the dim¬ 
ple shallower in the presence of constant the magnitude of of the quasilinear diffusion 
coefficient, D, must be larger by the factor Vi/{v‘f — 1). 

Nevertheless, when vt (the ion mach number) substantially exceeds 1 (the upper limit for 
ion-ion instability), the ion susceptibility contribution does not change the linear marginal 
stability condition by very much. It does not much enhance the required quasilinear diffusiv- 
ity nor the level of turbulence required to produce it, and it does not substantially raise the 
typical phase-space island width of the incoherent modes at which quasilinear stabilization 
occurs. 

Ion instability drive does not change the conclusion that incoherent quasilinear flattening 
of the electron dimple would occur at fluctuation levels that are too low to make signihcant 
non-linear modihcation to the ion distribution. The linear drive of the combined electron 
and ion distributions is brought to zero, if the ion velocities are signihcantly higher than the 
ion-ion stability threshold Cg, well before the quasilinear diffusivity of the ions is signihcant, 
and before entrainment of the ions into typical mode sizes. This conclusion is consistent with 
the code observation [T7] of a sustained initial period when the ion distribution evolution is 
quiescent and laminar, and the electron huctuations are predominantly localized to their 
phase-space separatrix. But it fails to explain the coherent structures that are observed 
to grow and entrain the ions even well before the ion beam velocities have slowed into the 
unstable regime. 

5 Electron Holes 

5.1 Hole Structural Relationships 

At the other end of the spectrum of treatments of non-linear effects and turbulence, far 
from the quasilinear diffusion approach, lies the phenomenon of phase-space “holes”. Con¬ 
ventionally such a hoie|2ni ED Ea Ea El ESI EH] refers to a localized coherent perturbation 
of the distribution function phase-space density that is self-binding via its self-consistent 
potential. A perturbation to a one-dimensional electron distribution function fe{y,v) can 
be self-binding if it traps electrons. To do so it must give rise to an electric potential that 
is positive in the hole. That requires the perturbation of the phase-space density, denoted 

^Incidentally, the electron-ion instability which slightly overlaps the ion-ion instability on a standard 
Stringer plot is suppressed by the flattening or hollowness of the electron distribution 
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/e to be negative: a deficit of electrons. (Likewise ion holes require a negative /*.) The 
parallel spatial coordinate is written y, for consistency with previous sections. Velocities v 
are unnormalized in this section. 

Electron holes in one dimension can exist at spatial scales from less than the Debye length, 
upward. The hole is essentially a Bernstein-Greene-Kruskal (BGK) mode|27]: a trapping 
structure that self-consistently satisfies the Vlasov-Poisson system of equations. There is 
substantial freedom in the form that such modes can take. Entropy arguments |23j support 
what is more often proposed as an ansatz [28] that the velocity dependence of fe in the 
trapped region is approximately parabolic (with positive curvature, negative temperature) 
leading to what is called a Maxwell-Boltzmann hole. They suggest that the most probable 
spatial extent of a shallow electron hole is approximately 4 times the plasma shielding length. 
However, the precise shape of the hole proves not to have a major effect on its properties, and 
modeling the hole as a rectangular box in y, v space, of constant depth, yields parameters that 
differ little from the Maxwell-Boltzmann hole|23]. It is therefore plausible to approximate 
the hole’s shape with simple model functions, and still expect to arrive at scalings that have 
reasonable quantitative validity. Moreover, deep holes with spatial extent much larger than 
the shielding length are possible. 

The self-consistent Poisson’s equation for an electron hole in one dimension of space may 
be written 



(30) 


where A is the plasma shielding length, and p is the charge-density of the hole. The shielding 
length would normally be thought of as the Debye length Xoe = \j€-QTe/ne^ but it can be 
generalized to account for ion shielding and for arbitrary distribution functions by regarding 
it as arising from the medium’s polarization term which is responsible for the dielectric 
susceptibility [23]. The real part of the linearized susceptibility for wave number k and phase 
velocity uj/k = Vp is then 




1 _ -p f 

k2\2 2^ p2 J dv Vp — 

species P 


(31) 


and this is the definition of A. fb denotes the unperturbed background distribution away 
from the hole; V denotes the principal value of the integral, and contributions like this from 
both electrons and ions are included. 

We consider a localized peaked hole potential structure that is the solution of Poisson’s 
equation (30). The potential energy then has a well, which Fig. [^illustrates schematically. 
For simplicity, and because it is the important case here, we take the hole to be stationary 
(corresponding to Up = 0), though moving structures can naturally be treated by a change of 
reference frame. It gives rise to phase-space orbits of electrons (along which fe is constant) 
that are the contours of constant kinetic plus potential energy, so they satisfy 


1 2 . 1 2 
-nieV + qcp = const. = -nieV^ 


(32) 


where we write q for the electron charge (it is negative), and Vb is the velocity at a distant 
unperturbed (“background”) position (yb) far from the hole, where the potential is 0 = 0 (so 
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Figure 5: Contours of constant / are contours of constant total energy in phase space. An 
attractive potential produces an island, which can become a hole. The parallel spatial extent 
of y here should be considered far less than the width of the wake shown in Fig. 


this potential is measured relative to the background plasma potential in the vicinity of the 
hole). Since the potential energy qcf) is negative, orbits that connect to yb have a minimum 
speed at any position given by 


Vsiy) = 


-2q(j){y) 


rrip 


(33) 


which is the boundary of the trapped-electron island in phase space. The charge density to 
be used in eq. (30) is then 

piy) = q J fiy,v)dv, (34) 

where / is the change of /. It is zero (/ = f{y,v) — f{yb,Vb) = 0, |n|> Us) on untrapped 
orbits , while on trapped orbits |n|< Vg, 


f{y, v) = f{y, v) - f{y, Vs) = f{y, v) - f{yb, o). 


(35) 


See Fig. The total electron density in the presence of the potential perturbation would 
be different at position y from its value at yb even if / were everywhere zero, because 


/ f{y^'^)dv 7 ^ / f{yb,v)dv. However, in Poisson’s equation (30), that difference is contained 
in the 0/A^ term, not in p. It is the linearized dielectric response of the plasma. In eq. (30) 
the p contains only the charge density attributable directly to /. 

Now we introduce a lumped-parameter model of the hole in which / has a characteristic 
magnitude at the center of the hole /o = /(O, 0), and the characteristic widths in velocity, v, 
and space, y, of the trapped region are An and Ay respectively. We dehne An so that the 
charge density at the hole’s spatial center {y = 0) is 


p = p(0) = g/oAn. 


(36) 


In Poisson’s equation the term d'^cp/dy'^ will be of magnitude approximately (p{0)/Ay‘^. But 
if the hole has large spatial extent. Ay 3> A, that term is negligible and we hnd as a quasi¬ 


lb 




























Velocity: v 


Figure 6: Schematic illustration of the distribution function f{v) in an electron hole. The 
reference distribution is uniform within the trapped region (|n|< n<j) and the difference / is 
what determines the binding charge density. 


neutral approximation to the hole 


A2 ^ eo 


-foAv. 

eo 


(37) 


(The quasi-neutral hole can be considered to be two “double-layers” that trap electrons 
between them. A small hole that is not quasi-neutral can be analysed [2^ |22] to hnd a 
comparable relationship between (j) and An, in which the proportionality coefficient depends 
upon Ay/X. So our conclusions are not qualitatively changed for holes of small Ay.) 

Hereafter, we refer to values at y = 0 and we drop the repetition of this fact in our 
notation. There is a proportionality between the two measures of the hole velocity width 
Vs and An. It requires knowledge of the velocity-shape of the hole to obtain its exact 
coefficient. For example, if / is uniform throughout the trapped region, then An = 2vs 
while if / is parabolic in the trapped region, then An = ^Vg, and if it is triangular An = Vg- 
Adopting this last alternative, for reasons that will become clear later, we have 


An = Vg 





and the relationship between the depth and velocity-width of the hole becomes 


-/ 


eofUeAv 

2g2A2 


rifeAn 

2a;2^A2 


AoeV n^Av 




MO) 


An 

Vte ' 


(38) 


(39) 


where Vte = ^j2Tf,/me and parameters such as Adb, ft and nt refer to the background 
plasma. Different assumptions about hole prohle shape will somewhat change the coefficient. 
But in general, when the shielding length is not too different from the Debye length, the 
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fractional hole depth f /ff, is roughly equal to the fractional velocity-width Av/vte for a 
quasi-neutral {Ay ^ A) hole. 

Since / cannot be negative there is a maximum hole depth and size obtained by setting 
the trapped distribution function (phase-space density) equal to zero i.e. —/ = /;,(0). Then 

Av _ 1 /A 

'^te \ ^De 



5.2 Hole Growth 

In view of the proportionality between the hole depth |/| and its velocity-width Vg, for a 
hole to grow in velocity width (and hence in potential) it must become deeper. When this 
happens in an effectively collisionless plasma, the shape of the hole, f{v), is determined 
not by maximizing entropy but by the constancy of / on orbits. The absolute value of / on 
trapped orbits {ft) is an invariant functional]] of the orbit’s action (/ vdx), provided the hole 
phase-space orbits remain closed. (Fine-scale mixing does not substantially change the mean 
/ on a phase-space orbit, and so does not escape this constraint. In recognizing it we abandon 
decisively the common presumption that the hole remains parabolic.) Therefore the only 
way for a hole to become deeper is that the external distribution function f{y,Vs) = f{yb, 0) 
increases. In Dupree’s analysis PU] of growth of moving {vp ^ 0) holes in an electron-ion 
instability, the way the external / increases is by the hole decelerating to lower \vp\ so 
that f{yb,Vp) increases (e.g. for a Maxwellian external distribution). In the present context, 
however, there is a different mechanism inducing hole growth. It is that the plasma is drifting 
in the x direction, perpendicular to the magnetic held. Consequently there is a convective 
time derivative of the external density: v^^', and this is indeed positive in the wake. In 
other words, the drift de-energization term in the Vlasov equation that has the effect of 
generating the dimple continues to operate if a hole is present, and is a cause of hole-growth. 
Or equivalently, perhaps conceptually simpler still, the hole experiences a background plasma 
of rising density because of drift. The importance of these remarks is that the growth of 
a hole is not suppressed quasilinearly by reaching sufficiently strong perturbation that the 
distribution function is flattened. The hole is coherent; and as long as that coherence is 
maintained, it continues to grow as f{yb, 0) grow^ 

The only way a hole stops growing short of maximal size, assuming there is insufficient 
turbulence to tear it apart, is for it to convect out of the spatial region where the external / 
is growing. Holes move mostly along the direction of the (unperturbed) phase space orbits. 
Therefore their parallel {y-) motion leaves them at approximately constant /. When a hole 
moves away from the peak of the wake’s potential prohle, it is therefore swept out of the 
wake, at approximately the electron parallel velocity, without any consequent growth. Once 
the hole reaches the unperturbed plasma outside the wake, no convective growth term is 
operating, and it will move away without further growth. 

Therefore, once a hole has formed with sufficient coherent integrity, the only condition 
for it to grow is that it stays inside the wake, which in general means it must remain near the 
peak of the wake’s potential energy curve (bottom of its electrostatic potential well) y = 0. 

here should be considered to be a distance large compared with the spatial extent of the hole but 
small compared with the width of the wake. 
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That hole position is unstable, so most holes will move from it, and then be convected out of 
the wake before they’ve grown large, because the electron orbit duration in the x-direction 
is rather small except when they are on axis. But a few may remain at the wake axis long 
enough to grow to near maximal size. It is those we now analyse. 

There is a linear relationship between / and Vg if A and the hole shape can be approx¬ 
imated as constant. Since the growing edge of a hole entrains additional phase-space area 
on which the distribution is equal to /^(O), the hole velocity prohle in this approximation 
is triangular. That was the basis for choosing the triangular prohle in the previous section. 
But we now do a self-consistent calculation that shows what the shape in velocity space of a 
growing quasi-neutral hole prohle actually is. This requires a treatment that is self-consistent 
and, for a deep hole, non-linear (i.e. avoiding the commonly used linearized plasma response 
explained in the prior section). It is most simply performed for a quasi-neutral hole by 
setting the net charge density to zero as follows. 

Consider a background plasma with Maxwellian electrons, 

fb{v)= exp(-u^), (41) 


where Vte = 2Te/me and u = v/vte- Write the normalized phase-space separatrix speed Ug at 
potential 0: = —q(f)/Te. This is positive because 0 is positive and q negative for electrons. 

Take the reference hat-top electron distribution to be constant within the trapped region, 
and constant along untrapped orbits as shown in Fig. 


„ . s _ Uf, J exp(M^ — M^) for u > Ug 
^ Vtey/^ ^ I ^ for M < Ms 


(42) 


This is the distribution that would arise if an electron-trapping potential hill arose slowly 
(compared with the electron bounce time) in a background distribution that was not varying 
with time. The density of this distribution is a function of the normalized potential, u^. It 
can readily be evaluated [20] as 


= J ffdv = rib 


^Ug r / s 

—= -I- e '’erlc(Ms) 


(43) 


The ions in the wake can be represented quite well [21 |T7| by two ion streams, each of 
narrow spread in speed. Near the wake axis they are equal and opposite. We take their Mach 
number outside the hole to be M = Vbjcg. Inside the hole the ion speeds are lower, because 
the hole repels ions, and by conservation of energy the mach number there is — 2ul. 

Hence, by conservation of flux, the ion density is 



M 


rib 


- 2ul 


(44) 


If the actual electron distribution in the trapped region is ft{v), different from the refer¬ 
ence // by /, then quasi-neutrality can be expressed as the cancellation of the charge arising 
from the background density in the perturbed potential, n/ — n*, and the hole charge density 
oc fdu. That is, 

rus 

Hf - Hi = nbH{Ug)/y/n = VtefbH{Ug) = 2vte / fb- ft du, (45) 

Jo 
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where in the trapped-region ff = fh is independent of u; fb refers to the background value 
/b(0) at M = 0; and the function H, for any given M, is essentially the normalized reference 
charge-density from both electron and ion distributions, 


H{us) = y/n 



-f e“'’erfc(-Us) 


M 

- 2ul 


(46) 


which must be cancelled by the hole charge-density. Figure]^ shows the form of H{u). 



Figure 7: The charge-density function H{u) for various values of ion mach number M. The 
solid line marked oo corresponds to an immobile ion background density independent of 
potential. 


Equation (45) determines the relationship between the hole velocity-width, u,,, and the 
changing background electron density expressed as the peak of its Maxwellian, fb- In the 
context of a growing hole, the actual trapped electron distribution ft remains invariant once 
formed. Therefore we can differentiate the equation with respect to the hole width, Ug, and 
it becomes 


d 

dug 





(47) 


No contribution arises to the differential from the fixed ft{u), and none comes from the limit 
because fb = ft sd u = Ug: the newly trapped electrons have phase-space density equal to 
the instantaneous background density. This equation can be written as a simple quadrature 


fdf_f dH 

J Y ~ J 2u-H' 


(48) 


The derivative H'{u) = dH/du is positive at moderate u, but reaches zero at a certain 
value u = Umax dependent on M. This is where the hole reaches its maximum possible size. 
At that size the rate of hole width-increase with respect to background density becomes 
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infinite: the hole blows up. Its growth can no longer be described by the quasi-neutral 
hole equilibrium equations. We call the corresponding value of the background distribution 
fb = fmax- And we hnd a set (at various M) of universal curves by numerical integration 
from Umax backwards toward zero. So that 

ln(///max) = f ^ ^ rr du. (49) 

'^max -tl 

Fig. shows the result. The signihcance of a curve for some given M is this. As a hole 



Figure 8: Trapped electron distribution function forms for various ion Mach numbers. These 
also give the hole growth curves arising in growing background density. 

moves from one value of / to a larger value, the hole size, Ug, grows by moving along 
the curve. As it does so, the trapped distribution function is built up within the hole by 
incrementally trapping additional phase-space. The shape of ft{u) is therefore also given 
precisely by this functional dependence. So these curves can be considered to represent 
the shape in velocity space [u) of a hole that grows from inhnitesimal size up to any hnite 
current size Ug- The boundary of the hole is Ug, at which ft = fb, and inside the hole 
ft{u) = fbexp[ln{f{u)/fmax) - ^^{fb/fmax)]- If the hole started at hnite size Uinu, with 
a trapped distribution inside Uinu different from the growing hole form, nothing changes 
except in the initial region of the hole u < Umu- It must begin from a trapped distribution 


that satishes the quasi-neutrality equation (45) but its initial shape will be determined by 


whatever mechanisms governed its formation. Thereafter, as it grows incrementally, driven 
by rising external density, it follows the curve, and the trapped distribution function is built 
up accordingly. 

No hole can grow stably beyond Umax, the place where the curves’ gradient becomes zero. 
That size is the maximum stable hole size. When it is reached, the hole blows up, and 
disruption of the ion streams will take place: their large free energy will be released through 
additional non-linear processes not described here. It is found that an excellent £t to the 
numerical values of Umax for the range of M shown is 

Umax = M/V 2 -O.Q. (50) 
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A hole that begins at a position with a certain valne of fb or eqnivalently backgronnd 
density rib, will grow as rib grows, provided it remains near the wake potential ridge, not 
convecting ont of the wake. It will reach the disrnptive size when the density has increased 
by a factor that can be read off the curves of Fig. For example, the M = 2.5 curve has 
ln(///maa;) = —1 at M = 0, SO a Small hole will reach disruptive size when the density has 
increased by a factor exp(l). Or again for M = 2, ln(///max) = —0.5 at u = 0.1, so a hole 
that starts with size Us = 0.1 will grow to disruption when n/riinit = exp(0.5). 

The wake axis experiences a large increase of background density {rib = '^oo exp(e0/Te)) 
as its potential subsides from the large negative values immediately behind the object. It 
is clear, therefore, that holes formed in this region that remain at the axis have sufficient 
density increase to reach disruptive size. This conclusion contrasts with the demonstration 
that quasilinear diffusion cannot reach a level of strong ion perturbation. 

These analytic conclusions are in accord with the observations in the simulations. Small 
electron holes form in the dimple (by mechanisms we don’t here calculate in detail). While 
they remain near the axis and are not convected out of the wake, they grow. Some eventually 
grow large enough to disrupt the ion streams. 


6 Summary 

The wake behind an object in a magnetized plasma with predominantly cross-held drift, 
and short Debye length, experiences one-dimensional electrostatic instabilities. The electron 
velocity distribution along the held acquires a depression we call the dimple, on orbits that 
spend a long time near the axial ridge of the potential energy structure of the wake. The 
driving term of this unstable dimple can equivalently be regarded as either de-energization 
of the electrons by drift perpendicular to B, down the potential energy ridge; or, as drift in 
an increasing background density, hlled in by parallel velocity less quickly for orbits with low 
parallel velocity near the ridge. The term may be approximated in the “Krook” collisional 
form in Boltzmann’s equation. 

The second viewpoint, of drift into increasing background density, provides a more trans¬ 
parent understanding of the resulting non-linear dynamics. The collisionless form of the dim¬ 
ple is immediately unstable to electrostatic waves near the wake axis, having phase velocities 
lying within the dimple velocity width of approximately {rn^jm^~^{^Cs = {m^jm^^l^Vte,j^2 
(approximately 7c<j for hydrogen plasmas). The waves will grow and, if incoherent, will hll in 
the electron distribution function dimple until it becomes marginally stable. The ion parallel 
velocity distribution consists of two streams attracted inward toward the wake axis. They 
contain a great deal of free energy. But, because the stream velocity spacing is large close to 
the object, the ion distribution is not itself linearly unstable until far downstream. Neverthe¬ 
less, it is observed in numerical simulations (and sometimes in space) that large-amplitude 
perturbations grow and substantially disrupt the ion streams, long before they have become 
linearly unstable. 

The electric helds associated with quasilinear velocity-space diffusivity sufficient to sta¬ 
bilize the electron distribution dimple are too small to cause substantial perturbation of the 
ions. The ion disruption therefore cannot be explained by linear or quasilinear instability 
growth. However, the driving term of electron instability is effective also in causing the 


21 


non-linear growth of electron holes. Holes formed away from the wake axis, or that move 
away from it, leave the wake at approximately the electron phase-space separatrix orbital 
speed before they can grow very much. Some holes, however, are formed at the wake axis 
and remain near it. When they do, the continuing background density enhancement causes 
them to grow to a maximum size beyond which they explode and disrupt the ion streams. 
We have calculated the form, the growth, and the maximum size of such electron holes, es¬ 
tablishing a (to our knowledge) new theoretical non-linear instability mechanism associated 
with cross-held drift into higher density regions. Electron holes grown by this mechanism 
provide the missing piece of the wake stability puzzle. Coherent hole growth is not sup¬ 
pressed quasilinearly, and explains how large, ion-disrupting, perturbations can occur before 
the ion streams become themselves linearly unstable. 

Undoubtedly, there are many other important phenomena in the moon wake, including 
Alfvenic processes that the present electrostatic treatment omits. But, because of their large 
linear growth rate, we consider the electrostatic instabilities to be primary. Several important 
details of the explanation remain to be investigated. We have not addressed the question of 
exactly how small electron holes form in the hrst place, nor have we quantitatively analysed 
their positional stability, which decides whether or not they remain at the wake axis and grow. 
Moreover, the present treatment, limited to parallel one-dimensional dynamics, omits oblique 
wave-vector perturbations, which might in some circumstances be important. Nevertheless, 
the qualitative agreement with phenomena observed in one-dimensional simulations provide 
strong evidence that the electron hole growth mechanism is the key explanation of these 
simulations at least. 
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